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ABSTRACT 



This paper investigates the application of the Kalman Filter and the 
General Exponential Smoothing techniques of forecasting. Both methods 
are derived and the similarities and differences between them are dis- 
cussed. The two techniques are then applied to the practical problem of 
predicting weekly losses suffered by the U. S. Marine Corps units in the 
I Corps Tactical Zone in the Republic of Vietnam. The mean absolute 
error of the prediction is used as the criterion for choosing the better 
of the two methods. Results are given for both techniques as well as 
for the method of linear regression. In general the Kalman Filter pro- 
vides the smallest mean absolute error for the three mathematical models; 
linear, growing sine with harmonics and frequency of sixteen, thirty-two, 
and fifty-two weeks, and a constant model. 



Library 

U.S. Naval Postgraduate School 
Monterey, California 93940 

TABLE OF CONTENTS 

I. INTRODUCTION - — — 7 

A. THE MODELS 8 

B. NOTATION — — 8 

II. THE KALMAN FILTER - 9 

A. THE KALMAN FILTER EQUATIONS 11 

B. REMARKS 13 

III. GENERAL EXPONENTIAL SMOOTHING — 14 

A. THE WEIGHTING VECTOR 15 

B. THE FORECAST BIAS — 17 

C. EXPONENTIAL SMOOTHING COMPARED WITH KALMAN FILTER 22 

IV. FORECASTING MARINE LOSSES 23 

A. FORECASTING WITH EXPONENTIAL SMOOTHING — — 24 

B. FORECASTING WITH THE KALMAN FILTER 42 

C. FORECASTING WITH LINEAR REGRESSION - 44 

V. CONCLUSIONS - 69 

COMPUTER PROGRAMS - 74 

BIBLIOGRAPHY - - 79 

INITIAL DISTRIBUTION LIST 80 

FORM DD 1473 81 



3 



LIST OF TABLES 



I. GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 

WITH HARMONICS, CYCLE LENGTH 16 WEEKS 27 

II. GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 

WITH HARMONICS, CYCLE LENGTH 32 WEEKS - 30 

III. GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 

WITH HARMONICS, CYCLE LENGTH 52 WEEKS - 33 

IV. GENERAL EXPONENTIAL SMOOTHING RESULTS, LINEAR MODEL 36 

V. GENERAL EXPONENTIAL SMOOTHING RESULTS, CONSTANT MODEL — 

VI. KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 

CYCLE LENGTH 16 WEEKS — - 45 

VII. KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 

CYCLE LENGTH 32 WEEKS 48 

VIII. KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 

CYCLE LENGTH 52 WEEKS - 51 

IX. KALMAN FILTER RESULTS, LINEAR MODEL - 54 

X. KALMAN FILTER RESULTS, CONSTANT MODEL 57 

XI. LINEAR REGRESSION RESULTS 60 



5 



I. INTRODUCTION 



One of the problems in any forecasting technique is what to do about 
additional data. That is, given a forecast based upon some data, how 
does one go about updating this forecast when additional data is obtained? 
Two such methods are discussed in some detail. These are the methods 
known as the Kalman Filter, named for R. E. Kalman, and General Exponen- 
tial Smoothing, a name coined by Robert G. Brown. 

Both the Kalman Filter and General Exponential Smoothing are derived 
and comparisons made for any similarities and differences. Both tech- 
niques are then applied to the problem of forecasting losses incurred by 
the U. S. Marine Corps operating in the I Corps Tactical Zone of the 
Republic of Vietnam. The actual losses are expressed in weeks and cover 
the period from December 24, 1965, to April 6, 1968. Losses include all 
Marine Corps personnel who die as a result of battle as well as non-battle 
causes. It also includes those personnel who must be evacuated from 
Vietnam because of wounds or injuries received from any cause. 

The basis for comparing the Kalman Filter and General Exponential 
Smoothing lies in their similarities. In each case, a mathematical model 
is assumed. The model consists of coefficients, which are unkown, and 
functions, called fitting functions, whose values are known. The problem 
becomes one of estimating the coefficients, based upon previous data, 
then using this estimate with the fitting functions to obtain a forecast. 
The forecast error, the difference between what actually occurs and what 
has been forecasted to occur, is weighed to obtain a new estimate of the 
coefficients. 
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A. THE MODELS 



The models used for forecasting U. S. Marine losses are a linear 
model , 

C(t) = a-j + a2t } 

a growing sinusoidal with harmonics, 

C(t) = (a-j + a£t) + (a^ + agt)sin~£- + 

(a* + a fi t)cos^i + a 7 sin^-^- + a Q cos^^- , 
v 4 6 ' w 7 w 8 w ’ 

and a constant model. 



C(t) 



In the second model, the forecasts were made for values of w of sixteen, 
thirty- two, and fifty- two weeks. 

The second model was chosen since there appeared to be some cycle in 
the number of losses incurred. It was felt that losses would increase 
as a result of increased enemy activity. Such activity was at a low ebb 
during the time of the monsoon rains. During the dry period, the tempo 
of enemy offensive actions increased for a time with a resultant increase 
in Marine Corps losses. 

B. NOTATION 

The notation used above reflects that used throughout this paper. 

ss /\ L% 

The symbols e(t) or C(t) represent the forecast of losses for the t Ln 
time period. The same symbols without the A above them are the observed 

j. L. 

value or the actual losses during the t Ln period. It is also convenient 
to adopt vector and matrix notation. Thus, an alphabetical symbol, such 
as a, will represent a column vector, while a T will denote a row vector. 
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In general, a vector with a subscript will represent the vector at a 
specific point in time. Matrices will be denoted by a capital alpha- 
betical character. In the event the matrix symbol is subscripted, it 
will denote the matrix obtained by taking the expected value of a col- 
umn vector, the first subscript, multiplied by the transpose of the vector 
denoted by the second subscript. A matrix with a single subscript is the 
expected value of the vector represented by the subscript times its trans- 
pose. For example, 

C gb = E(ab T ) , and 
C a = E(aa T ) . 

II. THE KALMAN FILTER 

Basic to every forecasting technique is the criterion which is used 
to determine the forecast. In the case of the Kalman Filter, which esti- 
mates the coefficients of the model, the criterion is to minimize the 

trace of the estimation error covariance matrix, C , where e is the 

e 

unobservable error vector between the true value of the model coeffi- 
cients and the current estimate of these coefficients. This estimate is 
the minimum mean square error estimator. 

A fundamental theorem of estimation, called the Gauss-Markov Theorem, 
is very important in any form of estimation. The theorem shall be stated 
but no proof is provided. 

Theorem: If o is the vector of the observed data at times t* and x 

is a linear function of q , x and are random variables with moment matrices 
C , C , and C , and if C"^ exists, then the linear minimum mean square 

X X0 

estimate x of x given t> is 



X 



and 



■ V 1 ! • 



c = c - c A cr 1 c T „ 

e x xe e xe 



Throughout this paper, the data will be assumed to be linearly 
related to the vector x (or a). Thus, 

0 = Bx + v, 



where v is the measurement error or noise vector and B is an m x n 
matrix whose rows are the vectors f^(t), the vector of fitting functions 
evaluated at time t. For example, if e(t) = x-| + x 2 t + v, the vector 
of fitting functions f^(t) = (l,t) and the vector of coefficients 
x — ( Xi , x 2 ) • 

Then, 



or 



T- 



X0 



X0 



E[x(Bx + v) ] 



C B' + C 

X XV 



and 



or 



c 0 = E[Bx, + v) (Bx, + y,) ], 
c e = BC x B T ♦ BC XV * cj/ ♦ C v . 



From the fundamental theorem, 



x = [C B T + C ][BC B T + BC + C T + C ] _1 e , 
L x xv JL x xv xv v J S ’ 



and C = C - [C B T + C ][BC B T + C T B T + BC + C ] _1 [C B T + C ] T . 
e x L x xv JL x xv xv v J L x xv J 



Normally the vectors x and v are assumed to be independent. This, 
with the assumption that E(v) = 0, implies that C is identically zero. 

A V 

With this assumption, the above equations reduce to 



x C x B T [BC x B T + C v ]o 
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and 



c = c - c b t [bc b t + c ] _1 [bc t ] t . 

e x x L x v x 

By the use of the matrix identities 

(2.1) AB T [C + BAB 1 ] -1 = [A -1 + B T C -1 B] _1 B T C -1 , 

and (2.2) [C _1 + B T A -1 B] -1 = C - CB T [BCB T + A] _1 BC, 

the following equations are obtained: 

x = [C" 1 + B T C~ 1 B]” 1 B T C ~ 1 0 

and C e = [C" 1 + B T C” 1 B] _1 . 

There still remains the matrix C which must be known before the 

X 

estimate x can be made. In general little information is available 
concerning x. In this case it is reasonable to assume that the diagonal 
elements of C are larger than any finite number. With this assumption, 

A 

(T 1 is considered to be 0. Then, 

X = [b T c~ 1 b] 1 b t c~ 1 e , 

and C = [B T C~ 1 B]* 1 . 

e L v J 

A. THE KALMAN FILTER EQUATIONS 

The foregoing has been the groundwork. Suppose that m observations, 
o m , have been taken and x m is the minimum mean square error estimate 
based upon the m observations. The matrix C™ is available. In addition, 
suppose that r additional observations are taken. The problem is to 

r +m 

determine the estimate x 

Since x m is available and is 

> = [ B ,nT ( C ,n ) - 1 B m ] - 1 b" ,T ( C m ) “ 1 



Also, 



"e 

m 



[B mT (C ra )-'B m ]-’ 



= B m x m + v m , 



and m+r _ R m+r m+r m+r . 

0 D X ' V 

m r 

The noise vector v is independent of the vector v . This implies that 
C™ +r is a diagonal matrix with elements C™ and C^. By partitioning the 

m-j- y* m r m*^ y* my* 

matrix B into components B and B and the vector e into e and e , 
the new estimate 

x m+r . [(c")- 1 ♦ B rT (cJ)-'B r r ,[(C e r ' xra + BrT ( C v>' , e r 3 

is obtained. Also, 

c m + r = [(c m ) -l + B rT (C^)- 1 B r ]- 1 . 

Then define K = c" +r B rT (cM"' . From equation (2.1), 



K = C^B rT [Cy + B r c”B rT r' 



From equation (2.2), 

„m+r 



C m - KB r C m 
e _ e 



Finally, the new estimate is 



"m+r 



= x m + K[o’ 



D r 'm- 
B x 



These last three equations are the Kalman Filter. [Liebelt 1967, 

pp. 165-166] Note that succeeding estimates x m+r are computed recursively 

from the prior estimate x m . 
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R. REMARKS 



Notice that with the assumption that C, is known, neither K nor (2 

v e 

depend unon the data e. These elements could be precomputed and re- 
tained for later use as needed. Also, it is not necessary to assume 
that the matrix C"^ is zero. If one has information about the distri- 
bution of x, the elements of C x can be comouted and used in determining 

C e* 

The vector K may be thought of as a weighting factor applied to the 
forecast error. Notice that the weight for any error is proportional 
to the inverse of the variance of the forecast. A large variance in 
the measurement or noise for any particular forecast will cause the 
error to be given less weight in determining the estimate of the coef- 
ficients to be used in computing the next forecast. 

r 

The dimension of the matrix C y is always r x r. If measurements 
or observations are taken at each succeeding time period, is a scalar 
quantity and is the variance of the measurement for that period. In the 
event that the matrix C v is diagonal with the diagonal elements equal to 
some constant raised to some power, the matrix C g will reach a steady 
state. That is , 

m + r .. m 
e e 

The power to which the constant is raised corresponds to the time period 
for which the observation is made. 
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III. GENERAL EXPONENTIAL SMOOTHING 



As is the case in the Kalman Filter method, the linear model 

x = Ae 

is assumed, but rather than estimating x in order to minimize the 
diagonal elements of some matrix, it is desired to minimize the sum of 
the weighted squared error, 

e w?[e(j) - e(j )] 2 . 
j-1 J 

If W is a diagonal matrix with elements (w^ w m ) , then 

£ w?[ e( i ) - e(i)] 2 = [W(e - e)] T W(e - e) . 

i-1 1 .... 

Also, it is assumed that 

/V /V 

0 = Bx , 

where the matrix B is defined the same as in the discussion of the Kalman 
Filter. 

Define the matrix F to be B^WW^B. Then 
[w(e - e)] T w( 0 - e) - [e T A T - e T WW T BF _1 ]F[e T A T - e T WW T BF -1 ] T 

- e T ww T BF" 1 [e T wwV + eVwo . 

The matrix F is positive definite. Therefore the choice of A^ which 
minimizes the above expression is 

a t = ww t bf _1 . 
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Fi nal ly , 



•'T T T -I 

x 1 = q'ww'bf 1 * 

Notice that this is identical with the Gauss-Markov estimate with 
F = B T C^B and WW T = C^ 1 

A. THE WEIGHTING VECTOR 

A fundamental hypothesis of Exponential Smoothing is that the 
vector of fitting functions f(t+l) must be a linear combination of f(t). 
The matrix of coefficients of this combination, which does not depend 
on time, is called the transition matrix, L. 



From this. 


f(t+l) 


= Lf (t) , 


and 


G(t+1 ) 


= f T (l )L T x(t-l ) , 


where 


x(t) = 


(L T ) t x(0) . 



The fittina functions are always computed at the number of time periods 
in the future for which the forecast is desired. This implies that the 
nresent period is always considered as the period t = 0. Then the 

A A 

problem becomes how to compute x(t) recursively from x(t-l) so that 

t _• ^ 2 

Z b [ 0 ( t - i ) - x(t)f(-i)] is minimized. 
i=0 ~ ~ 

Define the matrix 

F(T) = l b 1 f (-i )f T (-i ) , 
i=0 

and the matrix F to be the limit of F(t) as t becomes very larae. The 
matrix F will exist if each fitting function f^(t) is areater than the 
square root of b _t . 
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Let the data vector be 



g(T) = e b^ed-i )f(-i ) . 

1=0 

th T 

Then the i element of the vector g(T) is e b K e(T-k)f . (-k) . 

k-0 1 



Then 

and 



g(T) = e(T)f(0) + bL - 1 g ( T- 1 ) , 



F(T)x(T) = i b i f(-1)f T (-i)i;(T) . 
1=0 ~ 



But 



e(T-i) = f T (-Dx(T) . 



Hence , 



F(T)x(T) = g(T) . 



Substituting the steady state matrix F and combining the two equation 
for g(T), the expression 

x(T) = e(T)F _1 f(0) + bF" 1 L _1 Fx(T-l ) 

is obtained. [Brown 1963]. 

But bF'W = L T - F~ 1 f ( 0 ) f T ( 0 ) L T 



and 



f T (l) = f T (0)L T . 



Define the vector h to be F”'f(0). Then the expression for x(T) becomes 



x(T) = he(T) + [L T - hf T (l )]x(T-l ) , 

or x(T) = L T x(T- 1) + h[e(T) - e(T)] . 

This last equation is the Exponential Smoothing means of estimating the 
model coefficients. [Brown 1963, pp. 174-177]. The vector h is the weight 
of the forecast error applied in obtaining the new estimate x(T). Note 
that x(T) is computed in a recursive way from x(T-l). 



16 



B. THE FORECAST BIAS 

One may question the accuracy of the forecast obtained by means of 
Exponential Smoothing. It has been shown [Bessler and Zehna 1968] that 
for the model C(t) = a + e(t) , 

where E[C(t)] = a 

that the forecast is biased in the early stages by a factor of 1-b^, 
where b is the weighting scalar applied to past data. Since b is 
between zero and one, then E[C(t)] = a 

when b^ has vanished. However, this shows that even where the mathe- 
matical model is correct, there is an initial bias in the forecast. 

This may lead one to ponder what occurs when the assumed model is 
incorrect. It has been demonstrated [Brown 1963, p. 128] that in the 
case where the assumed model is 

E[C(t)] = a 

but the data is actually corresponding to the model 

E [C ( t ) ] = a ] + a 2 t , 

that E [C ( t ) ] = a-| + a 2 t - a 2 b(l-b)~^ , 

after b t vanishes. Thus, even after a sufficient amount of time has 
elapsed so that b^ = 0, the forecast is incorrect or biased by a constant 
times the slope of the line through the data. 

The next step was to determine if these results could be generalized. 
To do this, the Exponential Smoothing equation 

a(t) = iJa(t-l) + h[C(t) - f T (l)a(t-l)] 

was examined. This can be rewritten as 
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a(t) 



Then if 



= V [L T - hf T ( 1 )] k hC(t-k) + [L T - hf T (l )] t a(0) . 

k-0 " ~~ 

E[C(t)] = g T (t)a , 



E[C(t)] 



f T (l)V[L T -hf T (l)] k hg T (t-k)a + f T (l )[L T -hf T (l )] t E[a(0)]. 
k=0 



The above expression is true in general when the assumed model is 
C(t) = f T (t)a + e(t) 
and the data actually follows the model 

C(t) = g T (t)a + e(t) . 

In order to compare this expression with those obtained previously for 
the constant model , the values 

L = 1, 
f(t) = 1, 

and g(t) = 1 

are substituted to obtain 

t-1 . 

E [C( t ) ] = e [l-h] K h a . 

k=0 

This can be written as 

E[C(t)] = afl-b*] , 

where b is 1-h. Thus the result is the same as before. 

Next, it was desired to consider the situation where the data was 
modeled by E[C( t ) ] = a , but the data actually obeyed the model 
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E[C(t) = 8^ a^t. 



Then substituting the values 



f T d) - 1. 



and 



g T (t) = (1 ,t) , 
a T = (a 1 , a 2 ) 



into the expression for E[C(t)], 




a 9 ) T h . 



k=0 ' ‘ 

This becomes E [C ( t ) ] = (a^+a^) (1 -h*) - a 2 (l-h)h _1 (l-h t ) . 



which is identical with the prior result with b = 1-h . 

Notice that the expression for E [C( t ) ] has been used to determine if 
a bias exists when the correct mathematical model is assumed as well as 
to determine the amount of the bias when an incorrect model is assumed. 

It was desirable to determine if the forecast was biased when the assumed 
linear model was correct. The equation then becomes 



E[C ( t ) ] = (a 1 +a 2 t)f T (l)[I-L T +hf T (l)]‘ 1 [I-(L T -hf T (l)) t ]h 

- a 2 f T (l )[I-L T +hf T (l )] _1 [L T - hf T ( 1 )][I-l T +hf (1 )] _1 

[ I - ( L T - h f T ( 1 ) ) t_1 ]h- ( t-1 ) [I -L T +hf T ( 1 ) ] _1 [L T -hf T ( 1 ) ] t_ 1 h . 



When h*" vanishes, the expression becomes 



E[C(t)] = a^ + a 2 t - a 2 (l-h)h~\ 
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It became important to know if the vector [L^-hf T ( 1 )] t h vanished as 
t became large. If the vector fails to vanish, the forecast will remain 
biased. The vector will vanish if the eigenvalues of [L^=hf^(l )] are 
all between minus one and plus one. To demonstrate this, let 
A = L T -hf T (l), where A is n x n and nonsingular. Also, let the eigen- 
values of A be m.j, n^, • • . ,m n , where each m^ lies between minus one 
and plus one. Let v-| , . . . ,v be eigenvectors of A. These vectors 

form a basis for E n . 

Then Av^ = m^-j , 



Also, if x is in E n , then 

x = a 1 v 1 + ...+ a n y n . 



Premultiplying each side by A, 



Ax 



a,Av, ..... a n Av n . 



Substitute to obtain 



Ax = a,™,*, ..... a n m n v n . 

Repeating this process t times, 

A l x = a,m‘v, ..... a n m‘v n . 

Each itk = 0 for sufficiently large t. Therefore, 

aS< = 0. 

n + 

Since this is true for all x in E , A h = 0. 
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In the case of the linear model, the eigenvalues are equal to 
1 - (h-j + h^ ) 2” ^ . These eigenvalues will be between minus one and plus 
one if h-j + h 2 lies between zero and four. The h for the linear model 
is such that this is true, if b is between zero and one. To demonstrate 
this, recall that 

h = F _ 1 f ( 0) , 







00 

F = E 
1=0 


bVf-Df^-i 


and 




f T (t) = 


0, t) . 


Then 


F- 1 . 


7-b 2 


(l-b) 2 


« 




(l-b) 2 


(l-b) 3 b _1 



Then h T = (h^, h^) = [l-b 2 , (l-b) 2 ]. Hence, 

• h 1 + h 2 = 2(l-b) 2 , 

which is less than four and is greater than zero. 

Therefore, the vector [L^-hf^(l )]^h vanishes in the linear model for 
all values of b which are between zero and one. Then 

E[C(t)] - (a 1 +a 2 t)f T (l)[I-L T +hf T (l)]" 1 h 

- a 2 f T (l)[I-L T +hf T (l)]' 1 [L T -hf T (l)][I-L T +hf T (l)]" 1 h . 

Thus the forecast is unbiased if 

f T (l)[I - L T + hf T (l )] -1 h = 1 

and f T (l )[I-L T +hf T (l )] _1 [L T -hf T (l )][I-L T +hf T (l )] _1 h = 0 . 

For a specific example see Chapter V. 
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C. EXPONENTIAL SMOOTHING COMPARED WITH KALMAN FILTER 



The Exponential Smoothing equation is a special case of the Kalman 
Filter in the steady state. Recall that the matrix F was defined to be 
B T WW T B. If WW T is equivalent to the matrix C ” 1 used in the Kalman 
Filter equations, then F”^ is identical to the matrix C g derived in 
Chapter II. Thus, in the case of fitting functions consisting of poly- 
nomials, exponentials, and sinusoidals, which can be generated by a 
transition matrix, the Kalman Filter estimate becomes identical with 
that of Exponential Smoothing when C g reaches the steady state. However, 
since the Kalman Filter method does not require that the fitting func- 
tions be generated by means of a transition matrix, it may be applied to 
models other than those listed above. 

The primary difference between the Kalman Filter and Exponential 
Smoothing is the models for which they are valid. The Kalman Filter is 
valid both in a dynamic system as well as in the steady state. The 
Exponential Smoothing technique is strictly valid only in the steady 
state. 
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IV. FORECASTING MARINE LOSSES 



In order to see how the Kalman Filter and the Exponential Smoothing 
techniques compared, both methods, along with linear regression, were 
used to forecast U. S. Marine Corps losses in the Republic of Vietnam. 

The actual data used included one hundred and nineteen weeks from 
December 26, 1965, through April 6, 1968. The losses included those 
Marines who either were killed or were wounded and evacuated from 
Vietnam. 

Three models were investigated. These were a growing sinusoidal 
with harmonics with a cycle of sixteen, thirty-two, and fifty-two weeks, 

C(t) = a-|+a 2 t+(a2+a 5 t)sin(pt) + (a^+agt)cos (pt) 

+ a 7 sin(2pt) + a 8 cos(2pt) , 

a linear model , 

C(t) = a-j + a^t , 

and a constant model , 

C(t) = a . 

The objective was to minimize the mean absolute error where the absolute 

J. L « 

error for the t Ln forecast was the absolute value of C(t) - C(t). The 
mean error and the estimated variance of the errors were included for 
comparison. 

Both the Kalman Filter and Exponential Smoothing require an estimate 
of the coefficients a(0) in order to begin forecasting. Captain Paul 
William O'Brien, U. S. Marine Corps, had studied the problem of 
Exponential Smoothing in forecasting Marine monthly losses. [O'Brien 1968]. 
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Using a least squares line through the monthly losses from March 1965 
to December 1966 and the least squares Fourier fit through the deviations 

about the regression line, Captain O'Brien obtained initial values of the 
coefficients for the growing sinusoidal model. The coefficients of the 
least squares line through the data were used as the initial values in the 
linear model. These values were scaled by one-fourth and used as the ini- 
tial values for forecasting weekly losses. Thus the growing sinusoidal 
model was initialized with 

a T (0) = (-36.45, 15.675, 62.61, 31.3975, .6325, .6325, 

-35.1175, 54.76) , 

and the linear model with 

a T (0) = (-36.45, 15.675) . 

A. FORECASTING WITH EXPONENTIAL SMOOTHING 

In each model the value of the weighting factor by which the errors 
were discounted was chosen to be the solution to 

1 - b n = .25 , 

where n was the number of coefficients in the model. The weighting 
vector h used in the Exponential Smoothing equation was computed by 
determining the steady state matrix F [Brown 1963]. Then 

h = F" 1 f (0) . 

The transition matrix in each model was determined by inspection, 
rememberina trignometric identities for sin(A+B) and cos(A+B). 

The forecasting procedure involved several steps. These were: 

Step 1; Initialize the vector a(0) and compute f(l). 
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Step 2; Compute the forecast C(l) by 
CO) = f T (l)a(0). 

Step 3; Compute the forecast error 
e ( 1 ) = C(l) - CO). 

Step 4: Obtain a new estimate a(l) by 

a(l) = L T a(0) + he(l ) . 

Step 5; Compute the forecast C ( 2) by 
C(2) = f T (l)a(l). 

Step 6; Repeat Steps 3 through 5 for each successive forecast 
forecast with the appropriate time arguments. 

The first model investigated was the growing sinusoidal with 
harmonics. For this model the transition matrix is 







1 


0 


0 


0 


0 




0 


0 


0 










1 


1 


0 


0 


0 




0 


0 


0 










0 


0 


cos p 


sin p 


0 




0 


0 


0 






L 


= 


0 


0 


-sin p 


cos p 


0 




0 


0 


0 


1 


1 






0 


0 


cos p 


sin p 


cos 


p 


si n 


P o 


0 










0 


0 


-sin p 


cos p 


-si n 


p 


cos 


P o 


0 










0 


0 


0 


0 


0 




0 


cos 


2p sin 


2p 








0 


0 


0 


0 


0 




0 


-si n 


2p cos 


2p 




where 


p is 2 tt 


divided 


by the cycle length. 


Also 










f T 0) 


= 


(1, 


1 , sin 


p, cos p. 


sin 


P, 


COS 


p, sin 2p , 


COS 


2p) 



The values of h used were, for the sixteen week cycle 
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(0.064709, 0.0011416, 0.025939, 0.12583, 0.00046034, 
0.0022276, 0.020951, 0.059457), 



h 



T 



for a thirty- two week cycle 

h T = (0.072061,0.0012067,0.056675,0.12816,0.00096984, 

0.0021717, 0.042191, 0-049775), 

and for a fifty- two week cycle 

h T = (0.089418, 0.0013569, 0.10333, 0.13127, 0.0016449, 

0.0020314, 0.063796, 0.029308). 

The results obtained for this model for the three cycles are contained in 
Tables I through III. 

The second model was the linear model. The same procedure was followed 
as outlined above. The transition matrix was 



and 




0 

1 



f T (l) = (1, 1) , 
h T = (0.25, 0.01795) . 



The results from this model are contained in Table IV. 

The last model to which Exponential Smoothing was applied was the 
constant model. With this model 

L = 1, 

f T (l) = 1, 

and h = 0.25 . 
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WEEK 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



TABLE I 



GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 16 WEEKS 



LOSSES ESTIMATED LOSSES ERROR 



61 


47 


14 


64 


32 


32 


84 


30 


54 


53 


58 


-5 


87 


89 


-2 


116 


120 


-4 


119 


128 


-9 


108 


102 


6 


167 


55 


112 


405 


36 


369 


89 


107 


-18 


85 


107 


-22 


339 


139 


200 


208 


251 


-43 


80 


294 


-214 


116 


270 


-154 


141 


245 


-104 


133 


229 


-96 


125 


227 


-102 


158 


244 


-86 


164 


282 


-118 


82 


319 


-237 


203 


318 


-115 


139 


337 


-198 


133 


321 


-188 


198 


304 


-106 


151 


315 


-164 


140 


322 


-182 


175 


332 


-157 


359 


352 


7 


280 


405 


-125 


1253 


404 


849 


273 


627 


-354 


150 


491 


-341 


257 


348 


-91 


559 


307 


252 


258 


402 


-144 


290 


411 


-121 


382 


418 


-36 


547 


426 


121 


119 


451 


-332 


125 


353 


-228 


85 


317 


-232 


54 


341 


- 287 


343 


408 


-65 



2 / 



WEEK 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 



TABLE I (CONTINUED) 



LOSSES ESTIMATED LOSSES ERROR 



118 


556 


-438 


96 


580 


-484 


80 


541 


-461 


262 


467 


-205 


131 


447 


-316 


85 


408 


-323 


177 


395 


-218 


106 


434 


-328 


435 


444 


-9 


145 


512 


-367 


91 


439 


-348 


126 


338 


-212 


442 


274 


168 


86 


337 


-251 


121 


320 


-199 


153 


340 


-187 


327 


373 


-46 


191 


429 


-238 


220 


403 


-183 


336 


367 


-31 


280 


364 


-84 


237 


357 


-120 


202 


359 


-157 


300 


369 


-69 


273 


407 


-134 


389 


412 


-23 


646 


419 


227 


732 


460 


272 


359 


486 


-127 


400 


397 


3 


253 


362 


-109 


247 


336 


-89 


249 


350 


-101 


248 


373 


-125 


536 


376 


160 


249 


426 


- 177 


256 


363 


-107 


254 


324 


-70 


315 


333 


-18 


208 


406 


-198 


269 


464 


-195 


203 


525 


-322 


290 


521 


-231 


608 


495 


113 


357 


517 


-160 


260 


439 


-179 


189 


360 


-171 


166 


317 


-151 


262 


321 


-59 


216 


371 


-155 
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TABLE I (CONTINUED) 



WEEK LOSSES ESTIMATED LOSSES ERROR 



96 


224 


385 


-161 


97 


304 


366 


-62 


98 


493 


343 


150 


99 


234 


364 


-130 


100 


262 


323 


-61 


101 


466 


338 


128 


102 


367 


440 


-73 


103 


357 


500 


-143 


104 


414 


520 


-106 


105 


393 


505 


-112 


106 


243 


444 


-201 


107 


438 


337 


101 


108 


208 


319 


-111 


109 


404 


277 


127 


110 


948 


333 


615 


111 


735 


529 


206 


112 


1056 


585 


471 


113 


787 


656 


131 


114 


1339 


590 


749 


115 


453 


667 


-214 


116 


563 


499 


64 


117 


475 


449 


26 


118 


676 


442 


234 


119 


582 


517 


65 



Mean Absolute Error 


= 170 


Mean Error 


= -68.97 


Error Variance 


= 4.42 X 10 4 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



TABLE II 



GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 32 WEEKS 



LOSSES ESTIMATED LOSSES 



ERROR 



61 

64 

84 

53 
87 

116 

119 

108 

167 

405 

89 

85 
339 
208 

80 

116 

141 

133 

125 

158 

164 

82 

203 

139 
133 
198 
151 

140 
175 
359 
280 

1253 

273 

150 

257 
559 

258 
290 
382 
547 
119 
125 

85 

54 
343 



60 


1 


64 


, 0 


63 


21 


67 


-14 


65 


22 


77 


39 


101 


18 


125 


-17 


145 


22 


177 


228 


266 


- 177 


241 


-156 


210 


129 


246 


-38 


231 


-151 


176 


-60 


140 


1 


120 


13 


108 


17 


106 


52 


125 


39 


153 


-71 


162 


41 


212 


-73 


240 


-107 


264 


-66 


301 


-150 


316 


-176 


321 


-146 


331 


28 


386 


-106 


404 


849 


681 


-408 


619 


-469 


521 


-264 


466 


93 


505 


-247 


457 


-167 


433 


-51 


446 


101 


509 


-390 


444 


-319 


395 


-310 


346 


-292 


302 


41 



30 



49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 



TABLE II (CONTINUED) 



LOSSES ESTIMATED LOSSES 



ERROR 



118 


348 


-230 


96 


324 


-228 


80 


297 


-217 


262 


270 


-8 


131 


300 


-169 


85 


290 


-205 


177 


271 


-94 


106 


287 


-181 


435 


287 


148 


145 


385 


-240 


91 


389 


-298 


126 


381 


-255 


442 


388 


54 


86 


485 


-399 


121 


464 


-343 


153 


454 


-301 


327 


451 


-124 


191 


492 


-301 


220 


480 


-260 


336 


471 


-135 


280 


486 


-206 


237 


473 


-236 


202 


442 


-240 


300 


401 


-101 


273 


390 


-117 


389 


373 


16 


646 


391 


255 


732 


476 


256 


359 


564 


-205 


400 


524 


-124 


253 


496 


-243 


247 


429 


-182 


249 


374 


-125 


248 


334 


-86 


536 


310 


226 


249 


378 


-129 


256 


358 


-102 


254 


347 


-93 


315 


340 


-25 


208 


352 


-144 


269 


331 


-62 


203 


331 


-128 


290 


312 


-22 


608 


320 


288 


357 


415 


-58 


260 


420 


-160 


189 


394 


-205 


166 


353 


-187 


262 


317 


-55 


216 


321 


-105 
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TABLE II (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


96 


224 


320 


-96 


97 


304 


329 


-25 


98 


493 


366 


127 


99 


234 


451 


-217 


100 


262 


446 


-184 


101 


466 


445 


21 


102 


367 


495 


-128 


103 


357 


498 


-141 


104 


414 


489 


-75 


105 


393 


488 


-95 


106 


243 


474 


-231 


107 


438 


417 


21 


108 


208 


423 


-215 


109 


404 


366 


38 


110 


948 


379 


569 


111 


735 


544 


191 


112 


1056 


617 


439 


113 


787 


754 


33 


114 


1339 


776 


563 


115 


453 


926 


-473 


116 


563 


781 


-218 


117 


475 


678 


-203 


118 


676 


564 


112 


119 


582 


531 


51 



Mean Absolute Error = 161 

Mean Error = -74.73 

Error Variance = 3.92 X 10^ 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



TABLE III 



GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 52 WEEKS 

LOSSES ESTIMATED LOSSES ERROR 



61 


63 


-2 


64 


73 


-9 


84 


79 


5 


53 


87 


-34 


87 


83 


4 


116 


88 


28 


119 


102 


17 


108 


114 


-6 


167 


123 


44 


405 


147 


258 


89 


235 


-146 


85 


219 


-134 


339 


203 


136 


208 


260 


-52 


80 


269 


-189 


116 


238 


-122 


141 


220 


-79 


133 


211 


-78 


125 


200 


-75 


158 


187 


-29 


164 


185 


-21 


82 


185 


-103 


203 


161 


42 


139 


175 


-36 


133 


169 


-36 


198 


163 


35 


151 


177 


-26 


1 40 


176 


-36 


175 


174 


1 


359 


183 


176 


280 


244 


36 


1253 


273 


980 


273 


569 


-296 


150 


535 


-385 


257 


466 


-209 


559 


435 


124 


258 


491 


-233 


290 


451 


-161 


382 


426 


-44 


547 


432 


115 


119 


484 


-365 


125 


408 


-283 


85 


349 


-264 


54 


294 


-240 


343 


245 


98 
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L. 



46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 



LOSSES 



TABLE III (CONTINUED) 

ESTIMATED LOSSES 



ERROR 



118 


292 


96 


275 


80 


259 


262 


247 


131 


292 


85 


296 


177 


287 


106 


307 


435 


302 


145 


389 


91 


376 


126 


346 


442 


327 


86 


397 


121 


352 


153 


320 


327 


302 


191 


335 


220 


324 


336 


324 


280 


358 


237 


372 


202 


373 


300 


366 


273 


389 


389 


403 


646 


449 


732 


559 


359 


672 


400 


657 


253 


649 


247 


596 


249 


545 


248 


501 


536 


462 


249 


509 


256 


466 


254 


431 


315 


402 


208 


394 


269 


359 


203 


347 


290 


320 


608 


323 


357 


416 


260 


422 


189 


398 


166 


356 


262 


313 


216 


306 



-174 

-179 

-179 

15 

-161 
-211 
-110 
-201 
133 
-244 
-285 
-220 
115 
-311 
-231 
-167 
25 
-144 
-104 
12 
-78 
-135 
-171 
-66 
-116 
‘ -14 
197 
173 
-313 
-257 
-396 
-349 
-296 
-253 
74 

-260 

-210 

-177 

-87 

-186 

-90 

-144 

-30 

285 

-59 

-162 

-209 

-190 

-51 

-90 
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TABLE III (CONTINUED) 



WEEK LOSSES ESTIMATED LOSSES 



96 


224 


287 


97 


304 


275 


98 


493 


289 


99 


234 


355 


100 


262 


338 


101 


466 


331 


102 


367 


383 


103 


357 


398 


104 


414 


407 


105 


393 


429 


106 


243 


440 


107 


438 


406 


108 


208 


432 


109 


404 


389 


no 


948 


407 


111 


735 


576 


112 


1056 


653 


113 


787 


801 


114 


1339 


840 


115 


453 


1016 


116 


563 


901 


117 


475 


825 


118 


676 


730 


119 


582 


702 



Mean Absolute Error = 154 

Mean Error = -67.78 

= 3.94 X 10 4 



ERROR 

-63 

29 

204 

-121 

-76 

135 

-16 

-41 

7 

-36 

-197 

32 

-224 

15 

541 

159 

403 

-14 

499 

-563 

-338 

-350 

-54 

-120 



Error Variance 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 



TABLE IV 



GENERAL EXPONENTIAL SMOOTHING RESULTS 
LINEAR MODEL 

LOSSES ESTIMATED LOSSES ERROR 



61 


11 


50 


64 


42 


22 


84 


65 


19 


53 


88 


-35 


87 


94 


-7 


116 


109 


7 


119 


128 


-9 


108 


142 


-34 


167 


148 


19 


405 


170 


235 


89 


257 


-168 


85 


227 


-142 


339 


201 


138 


208 


258 


-50 


80 


260 


-180 


116 


222 


-106 


141 


203 


-62 


133 


196 


-63 


125 


187 


-62 


158 


177 


-19 


164 


179 


-15 


82 


182 


-100 


203 


159 


44 


139 


178 


-39 


133 


173 


-40 


198 


166 


32 


151 


181 


-30 


140 


177 


-37 


175 


171 


4 


359 


177 


182 


280 


236 


4 


1253 


257 


996 


273 


568 


-295 


150 


504 


-355 


257 


418 


-161 


559 


384 


175 


258 


449 


-191 


290 


406 


-116 


382 


382 


0 


547 


391 


156 


119 


448 


-329 


125 


360 


-235 


85 


295 


-210 


54 


233 


-179 


343 


177 


166 


118 


223 


-105 
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47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 



TABLE IV (CONTINUED) 
LOSSES ESTIMATED LOSSES 



ERROR 



96 


189 


80 


157 


262 


128 


131 


162 


85 


148 


177 


124 


106 


134 


435 


120 


145 


210 


91 


190 


126 


159 


442 


146 


86 


232 


121 


190 


153 


168 


327 


161 


191 


209 


220 


204 


336 


209 


280 


248 


237 


261 


202 


257 


300 


243 


273 


262 


389 


269 


646 


308 


732 


416 


359 


524 


400 


491 


253 


477 


247 


422 


249 


377 


248 


343 


536 


317 


249 


385 


256 


349 


254 


324 


315 


304 


208 


307 


269 


277 


203 


273 


290 


250 


608 


259 


357 


363 


260 


365 


189 


337 


166 


294 


262 


254 


216 


254 


224 


239 



-93 

-77 

134 

-31 

-63 

53 

-28 

315 
-65 
-99 
-33 
297 

-146 

-69 

-15 

166 

-18 

16 

127 

32 

-24 

-55 

57 

11 

120 

338 

316 
-165 

-91 

-224 

-175 

-128 

-95 

219 

-136 

-93 

-70 

11 

-99 

-8 

-70 

40 

349 

-6 

-105 

-148 

-128 

8 

-38 

-15 
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TABLE IV (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


97 


304 


231 


73 


98 


493 


249 


244 


99 


234 


321 


-87 


100 


262 


296 


-34 


101 


466 


286 


180 


102 


367 


340 


27 


103 


357 


351 


6 


104 


414 


357 


57 


105 


393 


377 


16 


106 


243 


387 


-144 


107 


438 


348 


90 


108 


208 


377 


-169 


109 


404 


330 


74 


110 


948 


353 


595 


111 


735 


536 


199 


112 


1056 


609 


447 


113 


787 


761 


26 


114 


1339 


794 


545 


115 


453 


984 


-531 


116 


563 


857 


-294 


117 


475 


793 


-318 


118 


676 


716 


-40 


119 


582 


718 


-136 



Mean Absolute Error = 128 



Mean Error 
Error Variance 



= - 2.21 
= 3.67 X 10 4 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 



ERROR 

61 

29 

37 

-8 

29 

49 

38 

16 

71 

290 

-102 

-79 

195 

14 

-117 

-52 

-14 

-18 

-22 

17 

19 

-68 

70 

-12 

-15 

54 

-7 

-16 

23 

201 

72 

1027 

-210 

-280 

-103 

225 

-133 

-67 

41 

196 

-281 

-205 

-194 

-176 

157 

107 



TABLE V 



GENERAL EXPONENTIAL SMOOTHING RESULTS 
CONSTANT MODEL 



LOSSES 


ESTIMATED LOSSES 


61 


0 


64 


35 


84 


47 


53 


61 


87 


58 


116 


67 


119 


81 


108 


92 


167 


96 


405 


115 


89 


191 


85 


164 


339 


144 


208 


194 


80 


197 


116 


168 


141 


155 


133 


151 


125 


147 


158 


141 


164 


145 


82 


150 


203 


133 


139 


151 


133 


148 


198 


144 


151 


158 


140 


156 


175 


152 


359 


158 


280 


208 


1253 


226 


273 


483 


150 


430 


257 


360 


559 


334 


258 


391 


290 


357 


382 


341 


547 


351 


119 


400 


125 


330 


85 


279 


54 


230 


343 


186 


1 18 


225 
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WEEK 



TABLE V (CONTINUED) 

LOSSES ESTIMATED LOSSES 



ERROR 



47 


96 


48 


80 


49 


262 


50 


131 


51 


85 


52 


177 


53 


106 


54 


435 


55 


145 


56 


91 


57 


126 


58 


442 


59 


86 


60 


121 


61 


153 


62 


327 


63 


191 


64 


220 


65 


336 


66 


280 


67 


237 


68 


202 


69 


300 


70 


273 


71 


389 


72 


646 


73 


732 


74 


359 


75 


400 


76 


253 


77 


247 


78 


249 


79 


248 


80 


536 


81 


249 


82 


256 


83 


254 


84 


315 


85 


208 


86 


269 


87 


203 


88 


290 


89 


608 


90 


357 


91 


260 


92 


189 


93 


166 


94 


262 


95 


216 


96 


224 



199 


-103 


173 


-93 


150 


112 


178 


-47 


166 


-81 


146 


31 


154 


-48 


142 


293 


215 


-70 


198 


-107 


171 


-45 


160 


282 


230 


-144 


194 


-73 


176 


-23 


170 


157 


209 


-18 


205 


15 


209 


127 


240 


40 


250 


-13 


247 


-45 


236 


64 


252 


21 


257 


132 


290 


356 


379 


353 


467 


-108 


440 


-40 


430 


- 177 


386 


-139 


351 


-102 


326 


-78 


306 


230 


364 


-115 


335 


-79 


315 


-61 


300 


15 


304 


-96 


280 


-11 


277 


-74 


259 


31 


266 


342 


352 


5 


353 


-93 


330 


-141 


295 


-129 


262 


0 


262 


-46 


251 


-27 
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TABLE V (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


97 


304 


244 


60 


98 


493 


259 


234 


99 


234 


318 


-84 


100 


262 


297 


-35 


101 


466 


288 


178 


102 


367 


332 


35 


103 


357 


341 


16 


104 


414 


345 


69 


105 


393 


362 


31 


106 


243 


370 


-127 


107 


438 


338 


100 


108 


208 


363 


-155 


109 


404 


324 


80 


no 


948 


344 


604 


111 


735 


495 


240 


112 


1056 


555 


501 


113 


787 


680 


107 


114 


1339 


707 


632 


115 


453 


865 


-412 


116 


563 


762 


-199 


117 


475 


712 


-237 


118 


676 


653 


23 


119 


582 


659 


-77 



Mean Absolute Error = 122 

Mean Error = 20.26 

Error Variance = 3.48 X 10^ 
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In this model the initial estimate of a(0) was taken to be zero. Also, 
since the initial forecasts are biased by 1-h*", the method of finite 
exponential smoothing [Bessler and Zehna 1968] was used. With this 
procedure, the weight given the forecast error in determining the new 
value of a(t) is h(l-h t )”^ rather than the constant h. This method 
increased the initial forecasts and as h* vanished, the forecasts 
became identical with those using the constant h throughout. These 
results are in Table V. 

B. FORECASTING WITH THE KALMAN FILTER 

The Kalman Filter requires, in addition to the initial value of 
a(0), the value of C”^ in order to forecast. In the models investigated, 
the matrix C v was not known. However, since the basic model was 

C = Ba + v 

and in forecasting one time period in advance, C v is a scalar, it was 
decided that could be estimated as the variance of the forecast error. 
Accordingly, the current estimate of the error variance was used as the 
value of C v in the Kalman Filter equations. 

The initial value of C g was obtained from the equation 

c e = [B T C~ 1 B ] - 1 . 

This required at least two forecasts in order to estimate C v - So that 
a better estimate could be made, C g was initialized after ten time 
periods. The forecasts for these ten periods were obtained by using 
the initialized value of a(0). 

The forecasting procedure was: 

Step 1; Initialize a( 10) = a(0). 
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Step 2; Obtain the forecasts for the first ten periods from 
C(t) = b T ( t)a(10) . 

Step 3; Compute the forecast errors, e(t), for these ten 
periods and the mean error, e. 

Step 4; Obtain an estimate of C v from 

1 10 9 

C = rg- z [e ( k ) - e]. 
v iu k=] 

Step 5; Initialize C g by 

C e (10) = [b t c” 1 b]” 1 . 

Step 6; Obtain new forecast C ( 1 1 ) by 
6(11) = b T (l l)a(10). 

Step 7; Compute forecast error e ( 1 1 ) and the new error mean, 
Step 8; Obtain new estimate of C y from 

1 11 9 

C v = yy z [e( k ) - eT . 

k=l 

Step 9; Compute weighting vector K by 

K = C e (10)b(ll)[C v + b T (ll)C e (10)b(ll)]" 1 . 

Step 10; Obtain new estimate a ( 1 1 ) from 
a ( 11 ) = a(l 0) + Ke(ll). 

Step 11; Compute new value of C e (11) from 

C e (ll) = C e (10) - Kb T (ll)C e (10). 



43 



Step 12; Repeat Steps 6 through 11 with the appropriate 

time arguments for each successive forecast. 

In general the value of C y will be different from C . An exception 

is in the constant model where they are identical. In this model the 

computation of K in Step 9 above involved the old estimate of the error 

variance as the value of C and C was the most recent estimate of the 

e v 

forecast error variance. 

C. FORECASTING WITH LINEAR REGRESSION 

The method of Least Squares was applied to the data using the 
linear model. The procedure employed was to fit the least squares line 
through the first twelve weekly loss figures. The coefficients of this 
line were used to obtain a forecast of losses for the thirteenth period. 
A new line was fitted to the thirteen observations and the new coeffi- 
cients used for the next forecast. This procedure was continued through 
each succeeding time period. 
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TABLE VI 





KALMAN FILTER 


RESULTS 






GROWING SINE WITH 


HARMONICS 






CYCLE LENGTH 1 


6 WEEKS 




WEEK 


LOSSES ESTIMATED LOSSES 


ERROR 


1 


61 


47 


14 


2 


64 


28 


36 


3 


84 


19 


65 


4 


53 


37 


16 


5 


87 


76 


11 


6 


116 


115 


1 


7 


119 


129 


-10 


8 


108 


107 


1 


9 


167 


58 


109 


10 


405 


10 


395 


11 


89 


-7 


96 


12 


85 


384 


-299 


13 


339 


-191 


530 


14 


208 


1146 


-938 


15 


80 


-269 


349 


16 


116 


-590 


706 


17 


141 


189 


-48 


18 


133 


539 


-406 


19 


125 


468 


-343 


20 


158 


285 


-127 


21 


164 


189 


-25 


22 


82 


100 


-18 


23 


203 


-61 


264 


24 


139 


80 


59 


25 


133 


91 


42 


26 


198 


131 


67 


27 


151 


231 


-80 


28 


140 


255 


-115 


29 


175 


250 


-75 


30 


359 


235 


124 


31 


280 


278 


2 


32 


1253 


242 


1011 


33 


273 


603 


-330 


34 


150 


441 


-291 


35 


257 


272 


-15 


36 


559 


221 


338 


37 


258 


322 


-64 


38 


290 


275 


15 


39 


382 


236 


146 


40 


547 


226 


321 


41 


119 


273 


-154 


42 


125 


191 


-66 


43 


85 


183 


-98 


44 


54 


223 


-169 


45 


343 


288 


55 



45 



46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 



TABLE VI (CONTINUED) 



LOSSES ESTIMATED LOSSES 



ERROR 



118 


435 


96 


453 


80 


411 


262 


339 


131 


340 


85 


330 


177 


342 


106 


389 


435 


387 


145 


422 


91 


321 


126 


197 


442 


118 


86 


178 


121 


180 


153 


220 


327 


261 


191 


306 


220 


275 


336 


234 


280 


237 


237 


256 


202 


293 


300 


332 


273 


374 


389 


366 


646 


339 


732 


331 


359 


329 


400 


268 


253 


258 


247 


253 


249 


265 


248 


272 


536 


260 


249 


293 


256 


265 


254 


269 


315 


312 


208 


393 


269 


448 


203 


484 


290 


464 


608 


419 


357 


404 


260 


346 


189 


301 


166 


283 


262 


287 


216 


309 


224 


301 



-317 

-357 

-331 

-77 

-209 

-245 

-165 

-283 

48 

-277 

-230 

-71 

324 

-92 

-59 

-67 

66 

-115 

-55 

102 

43 

-19 

-91 

-32 

-101 

23 

307 

401 

30 

132 

-5 

-6 

-16 

-24 

276 

-44 

-9 

-15 

3 

-185 

-179 

-281 

-174 

189 

-47 

-86 

-112 

-117 

-25 

-93 

-77 
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TABLE VI (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


97 


304 


273 


31 


98 


493 


257 


236 


99 


234 


290 


-56 


100 


262 


308 


-46 


101 


466 


363 


103 


102 


367 


459 


-92 


103 


357 


510 


-153 


104 


414 


513 


-99 


105 


393 


476 


-83 


106 


243 


406 


-163 


107 


438 


318 


120 


108 


208 


295 


-87 


109 


404 


274 


130 


110 


948 


310 


638 


111 


735 


409 


326 


112 


1056 


444 


612 


113 


787 


485 


302 


114 


1339 


463 


876 


115 


453 


503 


-50 


116 


563 


451 


112 


117 


475 


448 


27 


118 


676 


459 


217 


119 


582 


492 


90 



Mean Absolute Error = 166 

Mean Error = 11.33 

Error Variance = 6.28 X 10^ 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 



TABLE VII 






KALMAN FILTER RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 32 WEEKS 

LOSSES ESTIMATED LOSSES ERRCR 



61 


60 


1 


64 


63 


1 


84 


63 


21 


53 


61 


-8 


87 


62 


25 


116 


69 


47 


119 


82 


37 


108 


102 


6 


167 


127 


40 


405 


156 


249 


89 


184 


-95 


85 


-75 


160 


339 


179 


160 


208 


1155 


-947 


80 


617 


-537 


116 


969 


-853 


141 


-362 


503 


133 


162 


-29 


125 


246 


-121 


158 


183 


-25 


164 


173 


-9 


82 


117 


-35 


203 


-102 


305 


139 


215 


-76 


133 


153 


-20 


198 


165 


33 


151 


345 


-194 


140 


266 


-126 


175 


227 


-52 


359 


291 


68 


280 


608 


-328 


1253 


505 


748 


273 


1782 


-1509 


150 


1007 


-857 


257 


258 


-1 


559 


-12 


571 


258 


198 


60 


290 


32 


258 


382 


21 


361 


547 


144 


403 


119 


384 


-265 


125 


156 


-31 


85 


50 


35 


54 


-21 


75 


343 


-67 


410 



48 



46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 



ERROR 

16 

21 

33 

238 

42 

0 

104 

1 

311 

-95 

-187 

-182 

91 

-380 

-359 

-339 

-168 

-327 

-262 

-100 

-132 

-130 

-113 

37 

29 

162 

401 

408 

-52 

4 

-131 

-85 

-30 

19 

352 

43 

88 

112 

188 

74 

136 

47 

113 

386 

32 

-111 

-202 

-225 

-119 

-163 



TABLE VII (CONTINUED) 
LOSSES ESTIMATED LOSSES 



118 


102 


96 


75 


80 


47 


262 


24 


131 


89 


85 


85 


177 


73 


106 


105 


435 


124 


145 


240 


91 


278 


126 


308 


442 


351 


86 


466 


121 


480 


153 


492 


327 


495 


191 


518 


220 


482 


336 


436 


280 


412 


237 


367 


202 


315 


300 


263 


273 


244 


389 


227 


646 


245 


732 


324 


359 


411 


400 


396 


253 


384 


247 


332 


249 


279 


248 


229 


536 


184 


249 


206 


256 


168 


254 


142 


315 


127 


208 


134 


269 


133 


203 


156 


290 


177 


608 


222 


357 


325 


260 


371 


189 


391 


166 


391 


262 


381 


216 


379 
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TABLE VII (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


96 


224 


361 


-137 


97 


304 


340 


-36 


98 


493 


332 


161 


99 


234 


358 


-124 


100 


262 


334 


-72 


101 


466 


322 


144 


102 


367 


347 


20 


103 


357 


356 


1 


104 


414 


365 


49 


105 


393 


382 


11 


106 


243 


393 


-150 


107 


438 


374 


64 


108 


208 


384 


-176 


109 


404 


351 


53 


110 


948 


348 


600 


111 


735 


419 


316 


112 


1056 


444 


612 


113 


787 


500 


287 


114 


1339 


504 


835 


115 


453 


565 


-112 


116 


563 


504 


59 


117 


475 


460 


15 


118 


676 


414 


262 


119 


582 


404 


178 



Mean Absolute Error 


= 189.5 


Mean Error 


= 7,64 


Error Variance 


= 9,05 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



ERROR 

-2 

-10 

2 

-41 

132 

195 

105 

-92 

-105 

48 

-649 

-465 

249 

96 

158 

438 

516 

524 

694 

2155 

-8030 

-7363 

-5382 

-3635 

-2341 

-1515 

-1184 

-918 

-675 

-347 

-430 

656 

-941 

-649 

-141 

367 

16 

225 

429 

637 

194 

420 

622 

1043 

3494 



TABLE VIII 



KALMAN FILTER RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 52 WEEKS 

LOSSES ESTIMATED LOSSES 



61 


63 


64 


74 


84 


82 


53 


94 


87 


-45 


116 


-79 


119 


14 


108 


200 


167 


272 


405 


357 


89 


738 


85 


550 


339 


90 


208 


112 


80 


-78 


116 


-322 


141 


-375 


133 


-391 


125 


-569 


158 


-1997 


164 


8194 


82 


7445 


203 


5585 


139 


3774 


133 


2474 


198 


1713 


151 


1335 


140 


1058 


175 


850 


359 


706 


280 


710 


1253 


597 


273 


1214 


150 


799 


257 


398 


559 


192 


258 


242 


290 


65 


382 


-47 


547 


-90 


119 


-75 


125 


-295 


85 


-537 


54 


-989 


343 


-3151 
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TABLE VIII (CONTINUED) 



WEEK LOSSES ESTIMATED LOSSES 



46 


118 


4022 


47 


96 


1478 


48 


80 


960 


49 


262 


755 


50 


131 


803 


51 


85 


751 


52 


177 


698 


53 


106 


727 


54 


435 


708 


55 


145 


877 


56 


91 


836 


57 


126 


777 


58 


442 


747 


59 


86 


865 


60 


121 


788 


61 


153 


740 


62 


327 


710 


63 


191 


747 


64 


220 


720 


65 


336 


703 


66 


280 


720 


67 


237 


710 


68 


202 


684 


69 


300 


649 


70 


273 


641 


71 


389 


623 


72 


646 


631 


73 


732 


689 


74 


359 


750 


75 


400 


724 


76 


253 


707 


77 


247 


663 


78 


249 


624 


79 


248 


590 


80 


536 


561 


81 


249 


583 


82 


256 


559 


83 


254 


540 


84 


315 


525 


85 


208 


522 


86 


269 


504 


87 


203 


498 


88 


290 


481 


89 


608 


479 


90 


357 


527 


91 


260 


527 


92 


189 


509 


93 


166 


478 


94 


262 


443 


95 


216 


423 



ERROR 



-3904 
-1382 
-880 
-493 
-672 
-666 
-521 
-621 
-273 
- 732 
-745 
-651 
-305 
-779 
-667 
-587 
-383 
-556 
-500 
-367 
-440 
-473 
-482 
-349 
-368 
-234 
15 
43 
-391 
-324 
-454 
-416 
-375 
-342 
-25 
-334 
-303 
-286 
-210 
-314 
-235 
-295 
-191 
129 
-170 
-267 
-320 
-312 
-181 
-207 
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TABLE VIII (CONTINUED) 



WEEK LOSSES ESTIMATED LOSSES ERROR 



96 


224 


395 


-171 


97 


304 


368 


-64 


98 


493 


351 


142 


99 


234 


359 


-125 


100 


262 


328 


-66 


101 


466 


302 


164 


102 


367 


303 


64 


103 


357 


291 


66 


104 


414 


280 


134 


105 


393 


278 


115 


106 


243 


276 


-33 


107 


438 


257 


181 


108 


208 


270 


-62 


109 


404 


255 


149 


no 


948 


271 


677 


111 


735 


358 


377 


112 


1056 


414 


642 


113 


787 


506 


281 


114 


1339 


560 


779 


115 


453 


675 


-222 


116 


563 


677 


-114 


117 


475 


690 


-215 


118 


676 


689 


-13 


119 


582 


710 


-128 



Mean Absolute Error = 


651 


Mean Error 


-359.16 


Error Variance 


17.47 X 10 5 
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TABLE IX 



KALMAN FILTER RESULTS 
LINEAR MODEL 



WEEK LOSSES ESTIMATED LOSSES ERROR 



I 


61 


2 


64 


3 


84 


4 


53 


5 


87 


6 


116 


7 


119 


8 


108 


9 


167 


10 


405 


11 


89 


12 


85 


13 


339 


14 


208 


15 


80 


16 


116 


17 


141 


18 


133 


19 


125 


20 


158 


21 


164 


22 


82 


23 


203 


24 


139 


25 


133 


26 


198 


27 


151 


28 


140 


29 


175 


30 


359 


31 


280 


32 


1253 


33 


273 


34 


150 


35 


257 


36 


559 


37 


258 


38 


290 


39 


382 


40 


547 


41 


119 


42 


125 


43 


85 


44 


54 


45 


343 



-21 


82 


-6 


70 


10 


74 


26 


27 


41 


46 


57 


59 


73 


46 


88 


20 


104 


63 


120 


285 


135 


-46 


122 


-37 


127 


212 


185 


23 


209 


-129 


203 


-87 


203 


-62 


206 


-73 


206 


-81 


204 


-46 


208 


-44 


210 


-128 


200 


3 


209 


-70 


207 


-74 


203 


-5 


210 


-59 


209 


-69 


206 


-31 


208 


151 


234 


46 


247 


1006 


286 


-13 


294 


-144 


298 


-41 


305 


254 


324 


-66 


330 


-40 


337 


45 


348 


199 


366 


-247 


364 


-239 


361 


-276 


357 


-303 


351 


-8 
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TABLE IX (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


46 


118 


358 


-240 


47 


96 


355 


-259 


48 


80 


350 


-270 


49 


262 


344 


-82 


50 


131 


347 


-216 


51 


85 


343 


-258 


52 


177 


337 


-160 


53 


106 


335 


-229 


54 


435 


329 


106 


55 


145 


340 


-195 


56 


91 


336 


-245 


57 


126 


329 


-203 


58 


442 


324 


118 


59 


86 


335 


-249 


60 


121 


328 


-207 


61 


153 


322 


-169 


62 


327 


318 


9 


63 


191 


323 


-132 


64 


220 


321 


-101 


65 


336 


320 


16 


66 


280 


325 


-45 


67 


237 


327 


-90 


68 


202 


326 


-124 


69 


300 


324 


-24 


70 


273 


327 


-54 


71 


389 


328 


61 


72 


646 


334 


312 


73 


732 


353 


379 


74 


359 


374 


-15 


75 


400 


378 


22 


76 


253 


383 


-130 


77 


247 


382 


-135 


78 


249 


380 


-131 


79 


248 


378 


-130 


80 


536 


377 


159 


81 


249 


388 


-139 


82 


256 


386 


-130 


83 


254 


384 


-130 


84 


315 


382 


-67 


85 


208 


383 


-175 


86 


269 


380 


-111 


87 


203 


379 


-176 


88 


290 


375 


-85 


89 


608 


375 


233 


90 


357 


388 


-31 


91 


260 


390 


-130 


92 


189 


388 


-199 


93 


166 


384 


-218 


94 


262 


379 


-117 


95 


216 


377 


-i61 



TABLE IX (CONTINUED) 



WEEK 


LOSSES 




ESTIMATED LOSSES 


ERROR 


96 


224 






374 


-150 


97 


304 






372 


-68 


98 


493 






372 


121 


99 


234 






380 


-146 


100 


262 






377 


-115 


101 


466 






376 


90 


102 


367 






382 


-15 


103 


357 






385 


-28 


104 


414 






387 


27 


105 


393 






391 


2 


106 


243 






394 


-151 


107 


438 






391 


47 


108 


208 






396 


-188 


109 


404 






392 


12 


110 


948 






395 


553 


111 


735 






416 


319 


112 


1056 






429 


627 


113 


787 






450 


337 


114 


1339 






462 


877 


115 


453 






485 


-32 


116 


563 






488 


75 


117 


475 






494 


-19 


118 


676 






497 


179 


119 


582 






505 


77 


Mean 


Absolute Error = 


141 








Mean 


Error = 


-15.33 








Error Variance 


4.18 X 


10 4 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 



ERROR 

61 

64 

84 

52 

58 

39 

39 

5 

64 

242 

-75 

-32 

245 

72 

-128 

-92 

-8 

-9 

-12 

27 

7 

-76 

50 

-62 

-65 

33 

-47 

-56 

4 

184 

104 

993 

6 

-123 

-7 

301 

-13 

31 

95 

242 

-198 

-122 

-79 

-58 

266 



TABLE X 



KALMAN FILTER RESULTS 
CONSTANT MODEL 

LOSSES ESTIMATED LOSSES 



61 


0 


64 


0 


84 


0 


53 


1 


87 


29 


116 


77 


119 


80 


108 


103 


167 


103 


405 


163 


89 


164 


85 


117 


339 


94 


208 


136 


80 


208 


116 


208 


141 


149 


133 


142 


125 


137 


158 


131 


164 


157 


82 


158 


203 


153 


139 


201 


133 


198 


198 


165 


151 


198 


140 


196 


175 


171 


359 


175 


280 


176 


1253 


260 


273 


267 


150 


273 


257 


264 


559 


258 


258 


271 


290 


259 


382 


287 


547 


305 


119 


317 


125 


247 


85 


164 


54 


112 


343 


7 
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TABLE X (CONTINUED) 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


46 


1 18 


123 


-5 


47 


96 


118 


-22 


48 


80 


110 


-30 


49 


62 


96 


166 


50 


131 


141 


-10 


51 i 


85 


133 


-48 


52 


177 


120 


57 


53 


106 


176 


-70 


54 


435 


174 


261 


55 


145 


228 


-83 


56 


91 


166 


-75 


57 


126 


125 


1 


58 


442 


126 


316 


59 


86 


132 


-46 


60 


121 


90 


31 


61 


153 


120 


33 


62 


327 


140 


187 


63 


191 


141 


50 


64 


220 


191 


29 


65 


336 


194 


142 


66 


280 


198 


82 


67 


237 


270 


-33 


68 


202 


264 


-62 


69 


300 


242 


58 


70 


273 


299 


-26 


71 


389 


298 


91 


72 


646 


360 


286 


73 


732 


371 


361 


74 


359 


505 


-146 


75 


400 


401 


-1 


76 


253 


400 


-147 


77 


247 


391 


-144 


78 


249 


317 


-68 


79 


248 


266 


-18 


80 


536 


252 


284 


81 


249 


268 


-19 


82 


256 


250 


6 


83 


254 


255 


-1 


84 


315 


254 


61 


85 


208 


305 


-97 


86 


269 


304 


-35 


87 


203 


277 


-74 


88 


290 


254 


36 


89 


608 


290 


318 


90 


357 


290 


67 


91 


260 


356 


-96 


92 


189 


353 


-164 


93 


166 


302 


-136 


94 


262 


223 


39 


95 


216 


262 


-46 
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TABLE X (CONTINUED 



WEEK 


LOSSES 


ESTIMATED LOSSES 


ERROR 


96 


224 


262 


-38 


97 


304 


241 


63 


98 


493 


296 


197 


99 


234 


302 


-68 


100 


262 


254 


8 


101 


466 


261 


205 


102 


367 


267 


100 


103 


357 


355 


2 


104 


414 


356 


58 


105 


393 


400 


-7 


106 


243 


398 


-155 


107 


438 


391 


47 


108 


208 


438 


-230 


109 


404 


437 


-33 


110 


948 


406 


542 


111 


735 


415 


320 


112 


1056 


660 


396 


113 


787 


813 


-26 


114 


1339 


788 


551 


115 


453 


798 


-345 


116 


563 


581 


-18 


117 


475 


563 


-88 


118 


676 


548 


128 


119 


582 


637 


-55 



Mean Absolute Error = 109 

Mean Error = 41.33 

Error Variance = 9.23 X 10^ 
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WEEK 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

23 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 



TABLE XI 



LINEAR REGRESSION RESULTS 



LOSSES ESTIMATED LOSSES ERROR 



339 

208 

80 

116 

141 

133 

125 

158 

164 

82 

203 

139 
133 
198 
151 

140 
175 
359 
280 

1253 

273 

150 

257 
559 

258 
290 
382 
547 
119 

125 
85 
54 

343 

118 

96 

80 

262 

131 

85 

177 

106 

435 

145 

91 

126 



196 

252 

256 
224 
208 
201 
193 
184 

184 

185 
170 
179 

176 
172 
178 

177 
174 
177 
203 
217 

350 

351 
337 
336 
369 
366 
366 
376 
402 
383 
366 
347 
327 
334 
320 
306 

291 

292 
283 

271 
266 

257 

272 
265 
255 



143 

-44 

-176 

-108 

-67 

-68 

-68 

-26 

-20 

-103 

33 

-40 

-43 

26 

-27 

-37 

1 

182 

77 

1036 

-77 

-201 

-80 

223 

-111 

-76 

16 

171 

-283 

-258 

-281 

-293 

16 

-216 

-224 

-226 

-29 

-161 

-198 

-94 

-160 

178 

-127 

-174 

-129 
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WEEK 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 



TABLE XI (CONTINUED) 



LOSSES ESTIMATED LOSSES ERROR 



442 


248 


194 


86 


263 


- 177 


121 


253 


-132 


153 


246 


-93 


327 


241 


86 


191 


248 


-57 


220 


246 


-26 


336 


246 


90 


280 


253 


27 


237 


256 


-19 


202 


256 


-54 


300 


255 


45 


273 


259 


14 


389 


261 


128 


646 


270 


376 


732 


293 


439 


359 


319 


40 


400 


323 


77 


253 


330 


-77 


247 


329 


-82 


249 


327 


-78 


248 


326 


-78 


536 


324 


212 


249 


337 


-88 


256 


336 


-80 


254 


334 


-80 


315 


333 


-18 


208 


334 


-126 


269 


331 


-62 


203 


330 


-127 


290 


327 ' 


-37 


608 


327 


281 


357 


342 


15 


260 


345 


-85 


189 


343 


-154 


166 


339 


-173 


262 


334 


-72 


216 


333 


-117 


224 


330 


-106 


304 


327 


-23 


493 


328 


165 


234 


337 


-103 


262 


334 


-72 


466 


333 


133 


367 


340 


27 


357 


343 


14 
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TABLE XI (CONTINUED) 



WEEK 




LOSSES 




ESTIMATED LOSSES 


ERROR 


104 




414 




346 


68 


105 




393 




350 


43 


106 




243 




354 


-111 


107 




438 




352 


86 


108 




208 




357 


-149 


109 




404 




353 


51 


110 




948 




357 


591 


111 




735 




380 


355 


112 




1056 




395 


661 


113 




787 




421 


366 


114 




1339 




437 


902 


115 




453 




472 


-19 


116 




563 




474 


89 


117 




475 




480 


-5 


118 




676 




484 


192 


119 




582 




493 


89 


Mean 


Absolute 


Error = 


139 






Mean 


Error 


= 


9.9 







Error Variance = 4.53 X 10 



4 
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13 
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11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 
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WEEKLY MARINE LOSSES 



96 120 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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EXPONENTIAL SMOOTHING (CONSTANT) FORECAST WEEKS 
Figure 5 
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Figure 6 
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V. CONCLUSIONS 



In the case of the linear model, it was shown in Chapter III that 

E[C(t)] = (a 1 +a 2 t)f T (l)[I-L T +h f T (l)] _1 h 

- a 2 f T (l ) [I-L T +h f T (l )] _1 [L T -h f T (l )][I-L T +h f T (l)] _1 h 

for large t and when the actual data was such that 
E[C(t)] = a.| + a 2 t. 

Then, E[C(t)] = a-j + a 2 t + a 2< 

This showed that when the assumed model was correct, the forecast was 
biased by the slope of the data line. In the case of the problem of 
forecasting the Marine losses, the least squares line through the one 
hundred and nineteen points had slope 3.41. This suggested that Exponen- 
tial should have produced a forecast which is 3.41 more than the actual 
data, or an error of -3.41. The actual mean error was -2.21. Also, the 
Exponential Smoothing forecast would be above the linear regression fore- 
cast when the losses were increasing and below when they were decreasing. 
The results bear this out, although not around the time periods when the 
actual losses changed directions. For example, the losses increased 
during the time periods thirty-seven through forty. The Exponential 
Smoothing forecast with the linear model exceeded the linear regression 
forecasts during these periods. The losses then decreased during the 
periods forty-one through forty-four. The Exponential Smoothing forecast 
was above the linear regression forecast for the forty-first week but it 
then was below the linear regression forecasts through forty-four. 
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Error 

Standard 300 _ 
Deviation 
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ESI 6 
•4 


KF 

A 


LS 

A 




+ + 
ES32 ES52 


200“ 

ES* 


ES 

o 






100* 




KF 

<? 


Growing Sine Model 
Linear Model 
Constant Model 

i i 


1 


1 





-75 -48 -24 0 24 44 

Mean Error 



STANDARD DEVIATION OF ERRORS VS. MEAN ERROR 
Figure 7 
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The results obtained with the Kalman Filter and the growing sine model 
were generally comparable with the results with Exponential Smoothing. 

This is evidenced in Figure 7. The exception was the sine model with a 
cycle length of fifty-two weeks. The Kalman Filter forecast for this 
model fluctuated badly for the initial forecasts. However, the last 
seventy weeks of the data produced a mean absolute error of -321.71 and 
a mean error of -208.71. The last forty weeks produced a mean absolute 
error of 219 and a mean error -24.08. When these results are compared 
with the overall mean absolute error of 651 and the mean error of -359.16, 
it appears that the forecasts were improving as the effects of the initial 
forecasts were dampening out. Thus, the model whould not be totally dis- 
missed, especially since the Exponential Smoothing results with the cycle 
length of fifty- two weeks did not vary greatly from the results when a 
cycle length of sixteen and of thirty-two weeks was used. The Kalman 
Filter technique assumes that the forecast error variance is known. For 
the problem investigated, that variance was estimated, which was a variation 
on the Kalman Filter. If the true variance had been known, the results of 
all the models with which the Kalman Filter was used might have been im- 
proved. 

With the models used in forecasting Marine losses, the simpler the 
model, the better the results. Captain O'Brien found that the linear model 
produced better results than the more complicated growing sine model for 

the prediction of monthly Marine losses [O'Brien 1968]. The same was true 

o 

for the weekly forecasts. In addition, the even simpler constant model 
produced even better results. However, if the actual cycle period had been 
used rather than an estimate, the results of the Exponential Smoothing with 
the growinq sine model might have been greatly improved. The period length 
can be determined as described in R.G. Brown's book[Brown 1963, p. 66-75]. 
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In general, the estimated mean of the forecast errors was not zero. 
For example, the Kalman Filter with the constant model, on the average, 
forecasted forty-one fewer losses than actually occurred. With this 
information, a better policy would have been 

C'(t) = C(t) + u, 

where u is the mean forecast error, C'(t) is the estimated losses and 
C(t) is the forecasted losses from the method used. The variance of the 
forecast error will not be changed by this modification since C(t) would 
still be computed using C(t-l) - C(t-l). However, the mean error, where 
the error is 

e ' ( t ) = C(t) - C'(t), 

will be zero. 

Finally, the results indicated, as evidenced in Figure 7, that even 
when an estimate of the forecast error covariance matrix C v is used, the 
Kalman Filter results do not differ greatly from those of Exponential 
Smoothing. Since, except for the constant model, the Kalman Filter is 
more complicated than either the Exponential Smoothing or the Linear 
Regression methods, either of these latter methods should be used when 
C y is unknown or when no computer facilities are available. When C v is 
known and a computer is available, the Kalman Filter technique might 
provide better results than either of the other two methods. 

Under the criterion chosen, the minimum mean absolute error, the 
Kalman Filter method using the constant model was the best forecasting 
method. This method could be easily used for forecasting Marine losses 
in Vietnam. The calculations for this constant model are relatively 
simple and can be accomplished with a calculator without the aid of a 
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computer. The calculation of the estimate of the forecast error variance 

is the most complicated portion of the constant model. The mean absolute 

error and the variance of the forecast errors achieved with this model 

compare favorably with the data mean of 294.6 and variance of the data of 
4 

5.21 X 10 . Using this model would be more beneficial than merely using 
the data mean as the forecast of the losses for the next time period. 
Because of the difference between the two variances, one would be more 
certain of a minimum error with the Kalman Filter constant model. 
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GENERAL EXPONENTIAL SMOOTHING PROGRAM 



This program computes successive weekly forecasts for the model 

E[C( t) ] = a-j + a 2 t . 



c 

c 

c 



c 

c 

c 



DIMENSION A(2),FZERG(2),F0NE(2),H(3,2),B(3),XL(2,2) 
*XEST ( 1 19) , X( 11?), TIME! 119) ,U( 2) ,E( 119) ,RELER< 119) , 
*Y(2),R(2) 

DATA B/O. 86603, C . 9A868 ,C . 97A68/ 

LM= 1 19 

CO 91 J= 1, 119 
91 T I ME ( J ) = J 
DC 3 J= 1 , 2 
CC 3 1=1,2 
3 XL(I,J)=1.C 
XL( 1,2 )=0.C 
MMM=2 



READ(5,7) ( (H ( I , J ) ,J=1 ,MMM ) , 1=1,3) 

7 FORMAT ( AF 13 . 5 ) 

PE AD ( 5 , 6 ) ( X ( 1 ), 1 = 1, 119) 

6 FORMAT ( 11F6.C ) 

F0NE(1)=1.C 
FCNE(2)=1.0 
DO 21 L 1 = 1 ,MMM 
Y(LI ) = 0.0 
R ( LI ) = C .0 
21 CONTINUE 
DC 8 1=1,3 
E SCM = C .0 
A ( 1 ) =- 1 A5 • E 
A(2)=62.7 

CALCULATE ESTIMATE AND ERROR OF ESTIMATE 

DC 10 I T=1 »LM 
S = C • c 

CC AC LK= 1 ,MMM 
S=S+A(LK)*FCNE(LK) 

AC CONTINUE 
XEST { IT )=S 

E ( IT ) = X ( I T )-XEST { IT) 

ESUM=ESUM4E( IT ) 

EMEAN=ESUM/ IT 

RELERC IT) = E( IT) /XEST { IT ) 

UPCATE COEFFICIENTS 

CALL GMPRD(A,XL»Y, 1,2,2) 

CC 15 NL= 1 » MMM 
R ( NL ) = Y ( NL J 
15 CONTINUE 

DC 20 J= 1 , MMM 

UCJ) = 0.0 

U( J) = H( I ,JME( IT) 

20 CCNTINUE 

DC 25 M= 1 , MMM 
ACM) = RCM ) + UCM) 

25 CCNTINUE 
1C CONTINUE 

WRITE(6,10C) 1,8(1) 

IOC FORMAT Cl' 1 ,23X,'CASE ',11,' WEIGHTING CONSTANT, B = 



» 
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*F 10.5,//// ) 

CALL XMOMEME,LM, XMEAN , S IGMA , DE V ) 

SUMAB = 0.0 
DO 202 NN* 1 , LM 
SUMAB=SUMAE*ABS(E(NN) ) 

202 CONTINUE 
ABERR=SUMAB/LM 
WRITE<6.102) 

102 FCRMAT(20X, 'LOSSES' ,10X,' ESTIMATED LOSSES', 10X, 

♦ • ERROR *,7X, 'RELATIVE ERROR' ,///) 

DC 35 MN= 1 , LM 

WRITE (6, 103) MN,X(MN),XEST(MN),E(MN)»RELER(MN) 

103 FORMAT <3X, 13, 13X,F7.2,1CX,F8.2,16X,F8.2,1AX,F6.3,/) 
35 CONTINUE 

WRITE(fc,2C5)ABERR 

205 FORMAT (///,10X,'MEAN ABSOLUTE ERROR = ’,FA.O) 

WRITE(6,992)XMEAN, SIGMA, DEV 
992 FORMAT (/. 1CX,'MEAN ERROR = • , F 8 . 2 , // , 1 OX , 

♦ •ERROR VARIANCE = ', E 1A .7 ,//, 10X ,' STANDARD DEVIATION 

♦ ' *F 10 . 5 ) 

WR ITE ( fc, 20 1 ) 

201 FORMAT (1' 1 ,38X, 'PLOT OF TIME VS LOSSES AND ESTIMATED 
♦LOSSES',///) 

CALL FLOTP(TIME,XEST,LM, 1) 

CALL FL0TP(TIME,X,LM»3) 

WR I TE ( 6 ,203 ) 

203 FORMAT (' 1' ,50X,'PLCT OF ERRORS VS . T IME ' , /// ) 

CALL FLOTP (TIME, E,LM,0) 

8 CONTINUE 

FCNE(2) = FONE(2 ) ♦ 1 • C 
177 CONTINUE 
ST CP 
ENC 



SLBRCLTINE XMOMEN ( X , L , A , B ,C ) 
DIMENSION X ( 1 1 9 ) 

A=0. 0 

8 = 0.0 

DC 10 1=1, L 
1C A = A+ X ( I I 
A = A/L 
DC 20 1=1, L 
20 B = B+ ( X ( l )-A)^2 
B = B/L 
C = SCRT(B) 

RETURN 

ENC 
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KALMAN FILTER PROGRAM 



This program computes successive weekly forecasts for the growing 
sinusoidal model with a cycle length of 32 weeks. 



REAL*8 C , C INV 

DIMENSION A(8),E(8) f XL(8,8),E(119),X(119),XEST(119) 
*,RELER(119),U(8) 2 CE(8,8),CK(8,8),C(8,8),CINV(8,8), 
*TIME( 119),EB(25,8),Y(8),H(3) 

INTOBS=15 
L= 1 19 

RE AD (5,1 ) (TIME (I ) ,X(I ), I=1,L) 

1 FORMAT (F5.C»58X,F6.0) 

N = 8 

P=(2.0*3.141593)/32.0 
212 R EAO (5,12) ( H( I), 1 = 1, 8) 

12 FORMAT (4F8 *5 ) 

P P=2 • 0*P 
B( 1 ) = 1 . 0 
B(2) = l .0 
B ( 3 ) = S IN ( P ) 

B ( 4 ) = COS( P ) 

B ( 5 ) = S IN ( P ) 

B(6)=B(4) 

B ( 7 ) = S IN( PP) 

B ( 8 ) =CGS ( P P ) 

A ( 1)=-145.8/4.C 
A(2)=62. 7/4.0 
A(3)=250.44/4.0 
A ( 4 )= 1 25 • 55/4 . C 
A ( 5 ) = 2 .53/4 .0 
A ( 6 ) =A ( 5 ) 

A ( 7 ) =- 14C.47/4.C 
A ( 8 ) =2 19.04/4 .0 
DC 10 1=1,8 
DO 10 J=l,8 
10 XL ( I , J )=0 *0 
XL(1,1)=1.C 
XL(2, 1 )= 1 • C 
XL(2,2)=l*0 
X L ( 3 , 3 ) =CG S ( P ) 

XL ( 3 , A ) = S IMP) 

X L ( A , 3 )=-XL ( 3 , A ) 

XL ( A , A )=XL ( 3 , 3 ) 

XL(5,3 ) = XL (3,3 ) 

XL ( 5 , <i ) = XL ( 3 , A ) 

XL ( 5 , 5 ) = XL (3,3) 

XL ( 5 ,6 ) = X L ( 3 » 4 ) 

XL(6,3)=-XL(3,4) 

XL(6,4)=XL(4,4) 

XL(6,5)=XL(6,3) 

XL ( 6 , 6 ) = XL ( 6 » 4 ) 

XL ( 7 , 7 )=CCS ( PP ) 

XL(7,8) = SIMPP) 

XL ( 8 , 7 )=- XL (7,8) 

XL(8,8)=XL(7,7 ) 

E SLM=0 .0 
AB ERR=C .0 
DC 19 J=l, INTCBS 
SUM=0 . 0 
00 7 1=1,8 
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non 



7 



8 



901 



4 

19 



17 

18 

3 

444 

2 

ICC 



15 



16 



304 

303 

3C5 

20 

21 



23 

22 



24 



SUM=SUM*A( I ) *0 ( I ) 

XE ST I J)=SUM 
E(J)=X(J)-XEST(J) 

RELER( J) = E( J)/XEST (J) 

E SllM=E SUM+E ( J ) 

EMEAN=ESUM/J 
ABERR=ABEPP+AeS(E ( J)) 

SUMO.O 

00 8 1=1. J 

SUM=SUM+(E (J)-EMEAN)**2 
CV=SUM/J 
DO 901 1 = 1 » 8 
BB(J»I)=B( I) 

CALL GHPRD(XL»E»Y .8.8.1) 

DC 4 M = 1 , 8 
Bm = Y(M) 

CONTINUE 
DC 18 1=1,8 
DO 18 K=1 » 8 
SUM=C . 0 

DC 17 J=1 , INTOBS 
SLM=SUM+BB ( J,K)*BB( J, I) 

C ( K , I )=SUM/C V 

CALL GAUSS3(8,C,C,CINV,KER,8> 

GO TO (2,3) ,KER 
WR I T E ( 6 , 444 I 

FORMAT (////, 50X, 'MATRIX SINGULAR 
STOP 

00 100 1 = 1,8 
DC 100 J= 1 » 8 
CE( I, J)=CINV( I ,J) 

INT0BS=INT0ES*1 
DO 14 IT= I NTCB S » 1 19 



SLM=C .0 
DC 15 1=1,8 
SLM=SLM+B ( I )*A( I ) 

XEST( IT)=SLM 

E ( IT ) = X( IT )-XEST< IT) 

RELER ( IT)=E( ITI/XEST ( IT) 

ABERR=ABERR+ABS(E( IT) ) 

ESOM=ESUM+E (IT ) 

EMEAN=ESUM / I T 

SUM=0 .0 

CC 16 1 = 1, IT 

SLM=SUM+(E(I)-EMEAN)**2 

CV=SUM/IT 

DC 303 1=1,8 

SUM=C • C 

DC 304 J= 1,8 

SUM=SUM+CE(I,J)*B(J) 

U( I) = SUM 
SUM=0 .0 



DC 305 1=1,8 
SLM=SUM+B ( I )*U( I ) 

CC 20 J= 1 , 8 
CK(J)=U(J)/(CV+SUM) 

UPDATE ThE COEFFICIENTS 

DC 21 1=1, E 
A ( I ) = A ( I ) *CK ( I)*E(IT) 

DC 22 1=1,8 

SUM=C . C 

DC 23 J= 1 , 8 

SUM = SliM + B(J)*CE(J, I ) 

U( I )=SUM 

CC 24 1=1,8 

DC 24 J= 1 , 8 

CKCE ( I , J )=CK ( I )*U( J) 

CE(I, J) = CE(I, J )-CKCE( I, J) 
CALL GMPRD(XL,E,Y,8,8,1) 
DO 13 J= 1 , 8 



/// ) 
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14 CONTINUE 
WRITE<6.25) 

25 FCPMAT( , 1 , ,20X,‘LCSSES',10X, 'ESTIMATED LOSSES', 10X, 
♦•RELATIVE ERROR',///) 

♦•ERROR • ,7X, 'RELATIVE ERROR' ,///) 

ABERR=ABERR/L 
OEV=SQRT (CV) 

DC 26 MN= 1 , L 

26 WRITE (6,27 ) MN , X ( MN ) , X EST ( MN ) , E ( MN ) , RE LER < MN ) 

2 7 FORMAT (3X, 13, 14X,F6.C , 12X ,F7.2 , 18X ,F7 .2 , 14X, F6.3,/) 
WRITE(6,28) ABERR,EMEAN,CV,DEV 

28 FORMAT (10X, 'MEAN ABSOLUTE ERROR * *,F8.2,//, 

♦ •MEAN ERROR = • , F 8 . 2 , // , 1 OX , • ERROR VARIANCE = »,E15.7, 
♦//,10X, 'ERROR STANCARD DEVIATION = 'E15.7) 

WR ITE ( 6 ,29 ) 

29 FORMAT (• 1* ,38X, 'PLCT CF TIME VS. LOSSES ANC ESTIMATED* 
♦' LOSSES', //) 

CALL FLOTP(TIME»XEST ,L,L) 

CALL PLOT P(TIME,X,L,3) 

WR ITE (6, 30 ) 

30 FORMAT (• 1* ,50X,'PLCT CF ERRORS VS. TIME',//) 

CALL FL0TP(TIME»E,L,0) 

214 STOP 
END 
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